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We compute the three-loop corrections to the potential of two heavy quarks. In particular we 
consider in this Letter the purely gluonic contribution which provides in combination with the 
fermion corrections of Ref. [1[ the complete answer at three loops. 
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The potential between two heavy quarks constitutes a 
fundamental quantity in Quantum Chromodynamics. It 
enters in a variety of physical processes like the thresh- 
old production of top quark pairs and the description of 
charm and bottom quark bound states. Furthermore, it 
is crucial for the understanding of fundamental quantities 
of QCD, such as confinement. (See Ref. [2J for a recent 
review.) 

The idea to describe a bound state of heavy coloured 
objects in analogy to the well-established hydrogen atom, 
goes back to the middle of the 1970s [3|. Shortly after- 
wards, about 30 years ago, one-loop radiative corrections 
have been evaluated in the works [j, Q • It took almost 20 
years until the next order became available [fUl] which, 
at that time, was a heroic enterprize. The two-loop cor- 
rections turned out to be numerically quite important 
which triggered several investigations to go beyond. End 
of last year the fermionic corrections to the three-loop 
static potential have been completed [l|, H, [l(| ■ In this 
Letter we report about the pure gluonic part which com- 
pletes the three-loop corrections to the static potential. 

We present our results for the static potential in mo- 
mentum space where it takes the form 



V(\q\) = 
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Here, C A = N c and C F = (JV| - 1)/{2N C ) are the eigen- 
values of the quadratic Casimir operators of the adjoint 
and fundamental representations of the SU(N C ) colour 
gauge group, respectively, and a s denotes the strong cou- 
pling in the MS scheme. The one- and two-loop coeffi- 
cients a\ [3, and 02 [fiHl, El are given in Eq. (4) of 
Ref. [l[ where also the higher order terms in e, necessary 
for the three-loop calculation, are presented. In Eq. (JTJ 
we identify the renormalization scale /x 2 and the momen- 
tum transfer q 2 . The complete dependence on a can 
easily be restored with the help of Eq. (2) of Ref. [l]. 

A new feature of the three- loop corrections to V(|g |) 
is the appearance of infrared divergences [l~2j which is 
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FIG. 1: Sample diagrams contributing to the static potential 
at tree-level, one-, two- and three-loop order. Solid and curly 
lines represent quarks and gluons, respectively. In the case of 
closed loops the quarks are massless; the external quarks are 
heavy and treated in the static limit. 



represented by the ln(/x 2 /g 2 ) term in Eq . fll . It has 
been evaluated for the first time in Refs. [13, [lj] (see also 
Ref. [13); m Eq. fflwe adopt the MS scheme which has 
been used in Ref. [l4|. Let us mention that the infrared 
divergence cancels in physical quantities after including 
the contribution where so-called ultrasoft gluons inter- 
act with the heavy quark anti-quark bound state. An 
explicit result can, e.g., be found in Ref. [14| where the 
cancellation has been demonstrated in order to arrive at 
the measurable energy levels of the heavy-quark system. 
We note in passing that higher order logarithmic contri- 
butions to the infrared behaviour of the static potential 
have been computed in Refs. fltj [l7j. 

Before presenting our results for 03 let us provide some 
technical details. We generate the four-point quark anti- 
quark amplitudes with the help of QGRAF [18]. Some sam- 
ple diagrams up to three-loop order are shown in Fig. [T] 
In a next step they are processed further with q2e and 
exp QjHIiJ where a mapping to the diagrams of Fig.[2]is 
achieved. The mapping to two-point functions is possi- 
ble since the only dimenionful quantity in our problem is 
given by the momentum transfer between the quark and 
the anti-quark. Although there is only one mass scale 
in our problem technical complications arise from the 
simultaneous presence of static lines (zigzag lines) and 
rclativistic propagators (solid lines) which significantly 
increases the complexity of the reduction to master in- 
tegrals. For this task we employ the program package 
FIRE 211 in order to achieve a reduction to about 100 ba- 
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FIG. 2: One-, two- and three-loop diagrams. The solid line 
stands for massless relativistic propagators and the zigzag line 
represents static propagators. 






FIG. 3: Three-loop master integrals where the 0(e) part is 
only known numerically. The label "— iO" indicates that in- 
stead of the static propagator 1 / (po + iO) there is the propa- 
gator l/(po — iO). 



sic integrals, so-called master integrals. The latter have 
to be evaluated in an expansion in e which we achieve 
with the help of the Mellin-Barnes method (see, e.g., 
Refs. [22l - [2"5l |). We managed to compute all the necessary 
coefficients of the e expansion of the master integrals an- 
alytically with the exception of terms of order e 1 of the 
three diagrams shown in Fig. [3] Results for the master 
integrals as well as more details on their evaluation will 
be published elsewhere. As a crucial tool providing very 
important numerical cross checks of the analytical results 
we applied the program FIESTA (2(| which is a convenient 
and efficient implementation of the sector decomposition 
algorithm. The colour factors of the individual Feynman 
diagrams have been computed with program color [27| . 

In our calculation we allowed for a general gauge pa- 
rameter £ in the gluon propagator. For individual dia- 
grams we observe the appearance of terms up to £ 6 . We 
have checked that the coefficients of the £™ (n = 1, . . . , 6) 
terms are zero. 

In order to present our results we decompose the three- 
loop coefficient in the form 
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where ni is the number of light quarks and the first three 
coefficients on the right-hand side have been presented in 
Ref. [l| (see Eq. (6)). Whereas for the fermionic contri- 
butions there are seven different colour structures in the 
case of 03° there are only two. Note that the result of 
all colour structures containing a factor Cf are generated 
by iterations of lower-order contributions and thus do not 
contribute to a§ . Diagrammatically such contributions 
are easily identified since the corresponding Feynman in- 
tegrals contain so-called pinch contributions of the form 



l/(po + zO) x l/(po — iO) where po is the zeroth component 
of a loop momentum. Our result for ai°^ reads 
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Similarly to the fermionic contribution new colour in- 
variants appear which can be traced back to Feynman 
diagrams as the third one in the second row of Fig. [T] 
Expressed in terms of N c one has d a p cd d A bcd / 'N a — 
(iV c 3 + 6iV c )/48. The coefficient of df cd d a A bcd has already 
been presented in Refs. [H,[29j], the coefficient of the C\ 
term is new. 

Let us now discuss the numerical effect of the three- 
loop contribution to the static potential. Inserting the 
numerical results for the coefficients aj in Eq. ([1) we ob- 
tain 



V(\q\) 



4irC F a s (\q\) 
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where /1 2 = q 2 has been adopted in order to suppress the 
infrared logarithm and the ellipses denote higher order 
terms in a s . The term "209" in the three- loop coefficient 
receives a large contribution ("211") from the C\ term 
whereas the new colour structure only contributes with 
a coefficient "—2". From Eq. ((4]) we observe at one-, 
two- and three-loop order a large screening of the non- 
fermionic contributions by the n\ terms which is most 
prominent in the case of for n; = 5. Here the difference 
between a 3 °^ and the fermionic contribution is one order 
smaller than the individual pieces. 

In Tab. U we show the numerical evaluation of the 
square bracket of Eq. (|4]) for the charm, bottom and top 
quark case, i.e. for ni = 3,4 and 5, adopting the appro- 
priate values of a s . For charm the three-loop corrections 
are almost as big as the one- and two-loop contributions 
whereas for bottom the three-loop contribution is already 
a factor of four smaller than the two-loop one. In the case 
of the top quark one observes a good convergence: the 
three-loop term is already a factor ten smaller than the 
two-loop counterpart. 

As already mentioned above, V(\q |) itself is not a phys- 
ical quantity. Hence let us consider the ground state en- 
ergy Ei of a heavy quarkonium system which has been 
evaluated to the third order in perturbation theory in 
Ref. [30j where the contribution from 03 has been kept 
unevaluated. We are now in the position to complete the 
numerical analysis. It is convenient to write the pertur- 
bative contribution to E\ in the form 



E\ x - 



E\ + SE 



(i) 



SE 



(2) 



SE 



(3) 



(5) 



3 



m 


a; 


1 loop 


2 loop 


3 loop 


3 


0.40 


0.2228 


0.2723 


0.1677 


4 


0.25 


0.1172 


0.08354 


0.02489 


5 


0.15 


0.05703 


0.02220 


0.002485 



TABLE I: Radiative corrections to the potential where 
the tree-level result is normalized to 1 (cf. Eq. Q). In the 
second column we also provide the numerical value of a s cor- 
responding to the soft scale where p ~ m q a s (m q being the 
heavy quark mass). 



with the Coulomb energy = —Cpalm q /4:. m q is 
the heavy quark mass and the superscript in brackets 
indicates the order in perturbation theory. Adopting for 
the renormalization scale the choice /is = CFOi s {^s) rr >-q 
we obtain 
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= a 3 8 E?[ 104.82 + 3.186 
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= a s s E? I 83.386+ 1.473 



15.297 ln(a s ) , 



15.297 ln(a s 



15.297 ln(a s ) J , 
(6) 



where the contribution from 03 has been marked sepa- 
rately. One observes that the numerical effect amounts 
between 1 and 4% of the non- logarithmic constant. 

Finally, it is interesting to compare our results with the 
predictions obtained on the basis of certain assumptions 
on the perturbative expansion. In Ref. [31| a Pade ap- 
proximation in the coupling constant has been performed 
whereas the findings of Ref. [32j are based on renormalon 
studies. For /4 3 they predict 313 and 292, respec- 
tively, which overshoots the exact result by 40 to 50%. 

More recently, a detailed comparison of the pertur- 
bative result with lattice simulations has been per- 
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formed 33 1 with the aim to extract After trans- 

forming the result of Ref. [33| to momentum space us- 
ing the formulae provided in their Appendix one obtains 
202 < 4 0) /4 3 < 337. Thus the lower limit of the 
tively big) interval covers the exact result. 

To conclude, in this Letter the three-loop corrections 
to the static potential have been completed by evaluating 
the gluonic contribution. Our main result can be found 
in Eq. ([3]) where the three-loop coefficients are given for 
general colour structure. Numerical sizeable corrections 
are observed for the non-fermionic contributions which 
are partly canceled by the fermionic corrections evaluated 
in Ref. Q. 

Let us stress that the static potential constitutes a fun- 
damental quantity of QCD. It represents a building block 
in many physical quantities like the determination of the 
bottom quark mass from the T(1S) bound state or the 
third-order correction to top quark threshold production 
cross section at a future electron positron linear collider 
which would result in the most precise value for the top 
quark mass. The static energy is also a crucial object 
when comparing perturbation theory and lattice simula- 
tions (see, e.g., Refs. [33-36]). We also want to mention 
the extraction of the strong coupling constant from lat- 
tice simulations where again the static potential and in 
particular 03 plays an important role pa , [38[ . 

Note added: While finishing this paper we be- 
came aware of the preprint C. Anzai, Y. Kiyo and 
Y. Sumino, "Static QCD potential at three-loop order," 
I arXiv:091 1.43351 [hep-ph] , where also a 3 has been com- 
puted. We agree with their Eq. (10), however, obtain a 
better precision. 
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